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Abstract 

We study 1-loop MHV amplitudes in = 4 super Yang-Mills theory and in = 8 
supergravity. For Yang-Mills we find that the simple form for the full amplitude presented 
by Del Duca, Dixon and Maltoni naturally leads to one that has physical residues on all 
compact contours. After expanding the simple form in terms of standard scalar integrals, we 
introduce redundancies under certain symmetry considerations to impose the color-kinematics 
duality of Bern, Carrasco and Johansson (BCJ). For five particles we directly find the results of 
Carrasco and Johansson as well as a new compact form for the supergravity amplitude. For six 
particles we find that all kinematic dual Jacobi identities arc encapsulated in a single functional 
equation relating the expansion coefficients. By the BCJ double-copy construction we obtain a 
formula for the corresponding J\f = 8 supergravity amplitude. Quite surprisingly, all physical 
information becomes independent of the expansion coefficients modulo the functional equation. 
In other words, there is no need to solve the functional equation at all. This is quite welcome 
as the functional equation we find, using our restricted set of redundancies, actually has no 
solutions. For this reason we call these results virtual color-kinematics duality. We end with 
speculations about the meaning of an interesting global vs. local feature of the functional 
equation and the situation at higher points. 



1 Introduction and Main Results 



Originally proposed by Bern, Carrasco and Johansson [1, 2], the remarkable discovery of color- 
kinematics duality provides a powerful tool for exploring gravity amplitudes by starting from Yang- 
Mills amplitudes and "squaring" the kinematic numerators (double-copy construction). This tech- 
nique has led to a new way of thinking about tree and loop amplitudes. At tree level the validity of 
this construction has been proved [■ >] , and the existence of duality-respecting kinematic numerators 
has been explicitly shown for all numbers of external particles [1, 4-6]. At loop level, arguments 
from unitary cuts and soft limit strongly suggest that the double-copy construction should hold 
as long as the duality is achieved [2, 3, 7]. This together with the validity of a duality-respecting 
formula has be confirmed in various cases for 4 points up to 4 loops and for 5 points up to 2 loops 
[2, 8-11], the situation at 4 points 5 loops is currently under exploration [ ], and the self-dual sector 
in Yang-Mills has also been studied [16]. A natural next step is to study 1-loop MHV amplitudes 
with > 6 as well as the explicit connection between the maximally supersymmetric Yang-Mills 
theory and gravity here. The purpose of this paper is to give some first steps in this direction. 

In the 90 's, Del Duca, Dixon and Maltoni showed that one- loop amplitudes in Yang-Mills can 
be written as ^ [i'(] 

(1-1) A^n^ = 9ym E ^i<x2...<x„4'H^i^2 . . . cr„) 

where faia^-on is a combination of structure constants jahc contracted in a necklace-like form 
and J^n {p\02 ■ ■ - On) is the so-called leading color partial amplitude of the more standard color- 
decomposition [Is]. 

In [19, 20], a formula for J^n'{a\a2 ■ ■ ■ Cn) in terms of special integrals which gives physical results 
on all compact contours respecting the color-ordering was presented. We write it as 

[1.2) A]^'[aia2 . . . cr„j = ^ Pn[cricr2 • • • o-„j 

(cricr2 • • • o-n) 

where (ctict2 . . . o"„) is a shorthand notation for the standard Parke- Taylor denominator [ ]. In the 
rest of the paper we will not write the (super) momentum conservation delta functions and the 
coupling constant in any amplitude in order not to clutter equations. 
In this paper, we show that a simple combination of (1-1) and (1.2) 



A 

gives a formula that matches physical residues on all compact contours. 



^Here we use a notation A^^'' (or A^^^^) to indicate the level of the amplitude, where L is the number of loops, and 
L = refers to tree-level amplitudes. The curly letter denotes the full amplitudes and the ordinary letter denotes 
the partial amplitudes. The tilde sign indicates that the formula is written without helicity factor or SUSY delta 
functions. Yang-Mills amplitudes are denoted by the letter A (or A) and gravity amplitudes by M. 
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We propose to use Vn((Ti(T2 ■ ■ ■ Cn) as a basis of integrals not only in Yang-Mills but also for 
gravity amplitudes. The intention is to seek for a clearer observation of the relations between the 
two theories. In brief, these are totally "planar" objects enjoying cyclic and reflection symmetry as 
well as simple residues on all contours, and we propose to call them MHV polygons. 

By expanding the MHV polygons onto scalar loop integrals and studying the color-kinematics 
duality, we manage to recover Carrasco and Johansson's result for 5-pt case in A/" = 8 supergravity 
["■■], and further transform it back to our formulation, which is as simple and compact as that in 
Yang-Mills theory 

(1-4) Ml' = > r— rn(cxia2(T3(T4a5). 

.e53^xz.) (^i^2a3a4a5) eia^a,asa,) 



Here [cri(j20"3cr4cr5] denotes the Parke- Taylor denominator of an MHV tree amplitude while e(o"icr20"30"4) 
is the standard fully anti-symmetric tensor contracted with four vectors. 

Then as we go on to 6-pt case, we start to observe a new phenomenon. The formula for the 
gravity amplitude can still be expressed nicely as 

(1.5) Mf'= p^^2^^2!^P,(a,a,.,a,a,a,), 

with the coefficients 7 required to satisfy a set of constraints, each of which relates two 7's that 
differ only by a transposition of two adjacent labels; e.g. for the standard ordering, we have 

7(123456) ^ 7(213456) 



(123456)e(123456) (213456)e(213456) 
^ ■ ' [12] [34] [56] 312334^45^56(314 + 5|6][3|4 + 5|6) ^ 

(12) (34) (56) e(123456)e(213456)e(3456) 

where e{abcdef) = {ab)[bc]{cd)[de]{ef)[fa] — [ab]{bc)[cd]{de)[ef]{fa). For convenience, we refer to 
these constraints as a single functional constraint, in the sense that it depends on a choice of per- 
mutation of the labels. Quite surprisingly, all physical information in A^g becomes 7 independent. 
More exphcitly, once 'P6(c"iO"20"3C"4C"5cr6) is expanded, on the physical contour, in terms of only 
scalar box integrals the corresponding coefficients become 7 independent after using the functional 
equation. This fact is actually quite welcome as the functional equation has no solutions! 

The functional equation (1.6) is found by first expanding the MHV polygons in the Yang-Mills 
amplitude in terms of scalar integrals. We then add some particular set of redundancies to impose 
the color-kinematics duality conditions. It turns out that all conditions boil down to the single 
functional equation (1.6). Using double-copy construction we obtain the corresponding gravity 
amplitude. Just as for five particles, it is possible to transform it back to our formulation leading 
to the formula presented above. 

The fact that the functional equation has no solutions implies that the numerators we found in 
Yang-Mills which in principle satisfy the color- kinematics duality do not actually exist. This should 
not come as a surprise as our ansatz for the redundancies was not the most general one. What is 
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surprising is that even with a "virtual" solution to the color-kinematics duality conditions one can 
still square the numerators to get the corresponding gravity amplitude. 

The paper is organized as follows. In Section 2 we introduce MHV polygons and summarize 
their general properties, after which we prove the formula for n-pt 1-loop MHV amplitudes in super 
Yang-Mills constructed with these polygons. Then in Section 4 we analyze the 5-pt case to give a 
taste of the relation between this new formulation and the color-kinematics duality in the context 
discussed by Carrasco and Johansson. Afterwards, in Section 5 we focus on the simplest case where 
new phenomena would arise, where we analyze the color structure, obtain the condition on which 
color-kinematics duality is satisfied, and explain the "inconsistency" that appears. In Section 6 we 
go on to test the corresponding conjectured supergravity amplitude, where we provide evidences 
from quadruple cuts. The details of this analysis as well as a short introduction to color-kinematics 
duality are summarized in the appendices. In the end, we list out some possible future explorations 
along this line. 



2 General Properties of MHV Polygons 
2.1 Definition 

As has been stated in the introduction, we use chiral tensor integrals that possess unit leading 
singularities to build 1-loop MHV amplitudes [19, 20], in order to bring all physically important 
information into the coefficients attached to them. Specifically within the context of our discussion, 
we call these objects n-pt MHV polygons (or MHV n-gons). They are symbolically denoted as 
Vn{criO'2 ■ ■ ■ (Tn), and diagrammatically represented by shaded polygons (with the number of vertices 
showing in the middle). They are functions of the configuration of external particles (both the 
ordering and the kinematics data), and can be defined in two equivalent ways. One is that an 
n-pt MHV polygon shares the same configuration of leading singularities with its corresponding 
1-loop MHV Yang-Mills partial amplitude, but with the values of all non-zero leading singularities 
normalized. The other definition, as can be seen in (1.2), is by the following equation 



The equivalence of the two definitions is due to the fact that for any quadruple cut on the 1-loop 
MHV Yang-Mills amplitude, the associated leading singularity (whenever it is non-zero) is always 
the corresponding tree-level amplitude [22]. 

By definition, these objects should be invariant under both cyclic permutation and refiection of 
the sequence of the particle labels 

(2.2) 'Pn{cria2 ■ ■ ■ cx„_i(j„) = 'P„((T2(T3 . . . cr„cri), 

(2.3) Vn{ 

That is to say, for n particles, the length of the set {Vn} is (n — l)!/2, and each MHV polygon is a 
planar object with a fixed ordering of particle labels. 
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2.2 Reduction Formulas 



It is always possible to expand MHV polygons as a linear combination of scalar pentagon integrals 
and scalar box integrals, and the general method to obtain this expansion is worked out in [20]. 
Here we will only list out the 5-pt and 6-pt cases, which are needed in subsequent analysis. 
For any MHV pentagon (e.g. 7^5(12345)), there is one unique reduction formula 

T>,(12346) = fH^gp/,(l|2|3|4|5) 

^U|45](51>|13].....3,^( 

^ e 1234 4v I I I 

where the summation in the second line is performed over cyclic permutations of the sequence 
(12345); /„ denotes the ordinary scalar loop integrals with n loop propagators, and the vertical bars 
in the arguments separate external legs into groups that connect to different vertices on the loop. 
We also use the notation Sah = {Pa + Pb)"^, and the e symbol is defined by 

^2 e{abcd) = ^ie^upaP^aPbPcPd 

= {ab)[bc]{cd)[da] — [ab]{bc)[cd]{da) , 

and so it is completely antisymmetric in its arguments. 

Redundancies begin to appear in the reduction of MHV hexagons (e.g. "^6(123456)). From the 
method as given in [_u], by fixing label 1 and 6, we can obtain one specific formula 

7^e(123456) = - J] '"'"'"^f + '''^' /5(12|3|4|5|6) 



e(3456) 

>^ (12)[23](34)[41](1|2 + 3|4)^ 
^ 7(12M) /4(1|23|4|56) 



In the above expression, each summation is performed over cyclic permutations of even steps. And 
we use the following abbreviations 

(2.7) {a\b + c\dy = {a\b + c\d][a\b + c\d) , 

(2.8) e{abcdef) = (ab) [be] (cd) [de] (e/) [fa] - [ab] (be) [cd] (de) [ef] (fa) . 

In (2.7) the square is just a notation to keep track of the dimension. We can observe that e{abcdef) is 
invariant under cyclic permutation of even steps, and acquires a minus sign under cyclic permutation 
of odd steps. 
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Note that, (2.6) is manifestly invariant only under cyclic permutations of the labels by even 
steps, and it consists of only half of all the pentagon and box integrals with the correct ordering. 
Although cyclic permutations of odd steps also leads to a valid formula (which is by itself also 
manifestly invariant under cyclic permutations of even steps), their equivalence is only guaranteed 
by non-trivial identities among the pentagon and box integrals. Let us denote these two reductions 
as 7/(123456) and 7/(234561) respectively. With the purpose of exploring color-kinematics duality 
later on in the context of loop integrals, we should exhaust the entire loop integral basis. And so 
we need to express the MHV hexagon as a linear combination of the two reductions 

(2.9) ^6(123456) = a(123456)7/(123456) + /3(123456)'H(234561), 
under the constraint 

(2.10) a(123456) + /3(123456) = 1. 

Here the argument in the parenthesis is pure labeling and identical up to cyclic permutations of 
even steps. If we regard them as functions of particles' kinematics data, the functions should respect 
this symmetry. But notice that by now there is no definition for e.g. q;(234561) and /3(234561), so 
that we can further identify 

(2.11) a(234561) = /3(123456), 

and require that as functions, q;(234561) and q;(123456) are related by permutations. Then both the 
a functional parameter and the reduction formula (2.9) manifestly enjoy the full cyclic invariance. 



3 1-Loop MHV Amplitudes in A/" = 4 Super Yang-Mills 

We first go on to show that given the MHV polygons defined in the previous section, the formula 
(1.3) 

matches correct physical residues on all contours. In order to do this, we need to use a special 
formula for the tree-level amplitude [2-^>], which can be easily proved by Britto-Cachazo-Feng-Witten 
(BCFW) method [24, 25] 

/q 1 \ T^'^) — /lo-2,0-3,...,a-„_2,o-n-l^ 

where we take the abbreviation 

(3.2) /i CT2,cr3,...,(T„_2,o'„_in /lo-2a2/a2cr3a3 ' ' ' /a„_2fn-l™' 
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and the permutation (Xj is taken over the label set {2, 3, . . . , n — 1}. This can be visualized as the 
summation of chain-like diagrams where we fix the two ends of the "chain" (in the given formula 
we fix 1 and n) and fully permute all the intermediate vertices. 

Then we go on to evaluate the factorization of 1-loop full amplitude in the quadruple cuts in 
A/" = 4 super Yang-Mills. Since in (3.1) each color factor is accompanied by a kinematic factor 
which has exactly the form of an MHV tree-level partial amplitude, we can start by evaluating the 
quadruple cuts on partial amplitude. As has been mentioned previously, in MHV super Yang-Mills 
amplitudes, whenever the factorization is non-trivial, it always has the form 

4 

(3.3) AW^l[A^^) = A^^).Bet{\J\), 

where Det(| J|) is exactly identical to the Jacobi determinant arising from cutting loop propagators 
and is to be exactly canceled by that factor under quadruple cut [22]. So what is left to be matched 
from the ansatz side is purely the Parke- Taylor form. 

Now switch to the full amplitude, where the non-trivial factorizations only come from two types 
of quadruple cuts. The virtue of the formula (3.1) is that, due to the freedom in picking up any 
two external particles as the two ends of the chain, if we always choose in the factorized amplitudes 
those "external particles" from the cut propagators, then as a result the color factors will just 
glue together to form the maximal loop (in the sense of color diagrams). And when combining 
the results on the corresponding kinematic factors that mimic partial amplitudes, we conclude the 
non-vanishing factorizations are always 

^„ ^ii^^ -i^etuju 2^ 2^ 2^ 2^ a ami ami ami am)' 

where the summation is performed over the group of external particles attached to each factorized 
amplitude respectively. This is exactly what we would get if the formula sums over all non-equivalent 
MHV polygons and the coefficient in front of each MHV polygon has the form as shown in (1.3). 
So we conclude that with the MHV polygons defined in (2.1) as the fundamental building block, 
the formula for 1-loop ra-pt MHV super Yang-Mills full amplitude is 

(3.5) Yl JT"Z V^{a,a,...a^). 



4 From Yang-Mills to Gravity: 5-pt 1-Loop Amplitudes 

In order to provide an example of how our formulation may work between Yang-Mills and gravity, 
in this section we focus on the 5-pt case, showing that our formula 
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naturally gives rise to color-kinematics duality in its expansion onto scalar loop integrals, and further 
obtaining the corresponding gravity amplitude as a compact expression also built purely upon MHV 
pentagons. 

In the reduction formula for the MHV pentagons (2.4), it is easy to observe that each receives 
a unique contribution from its corresponding "Ps, while the coefficient in front of each I 4 would 
receive contributions from two Ps's. To make the structure that appears in later discussion more 
apparent, we define 

(4.2) Qia.a^a.a.a,) = h^^^, 

e(cri(T20-30-4) 

where similar to the abbreviation for angle brackets, we have 

(4.3) [(Tl(T2 ■ ■ ■ Cr„] = [criO-2] [cr20-3] ■ ■ ■ [(T„(Ti] . 

The function Q is totally symmetric in any cyclic permutation of the labels and acquires a minus 
sign under refiection ^ (This is desirable since {(yia2(r3(y iO'^)Q{aia20'3a ^a^) is the coefficient of the 
pentagon integral, and we would prefer it to have the same symmetry with the corresponding basis 
element). Then the coefficient of one specific /4, e.g. /4(1|2|3|45) is 

[12] [13] [23] [45] [12] [13] [23] [54] 

^^^^ (45)e(1234) ^™ (54)e(1235) 

= — f alb fb2cfc3df deaf e45 [Q(12345) - Q(12354)] . 

"545 

So in the expansion the formula is again purely in terms of another functional coefficient Q 

■A.'^^ = /<7i(T2<73,f74cr5 [Q{(^l<^2<y3<74<^5) - Q(o-lCr20"30"5C^4)] ^ 
(4.5) ^ ^-405 

02 03 0-40-5 

Q(cri(J2Cr3Cr4(T5)/5(cri|o-2|o-3|cr4|(T5), 

where we use the abbreviation 
(4.6) 

(T20-3, 0-40-5 f aaibfbcr2cf ccrgd f deaf €(74(75 j 

and the summations are performed over all nonequivalent permutations of the labels respectively. 
Very nicely, at this point we may observe that this is exactly Carrasco and Johansson's result 
obtained in [ ■] from an ansatz that respects color-kinematics duality where the basis are {/4/s, I^}. 

A more interesting implication is that the MHV pentagons in the new formulation naturally 
encodes the color-kinematics duality in an implicit way. In fact, we can modify the expansion of V5 

(2.4) into a different form, e.g. 

(4.7) ^5(12345) = (12345)|g(12345)/5(l|2|3|4|5) + ^^^^^^^^ ~ ^^^^^^^^ /^(i|2|3|45) 



•^It is interesting that this Q also played an important role in the analysis of [s] 
■^For a quick revie-w of the color-kinematics duality, please refer to Appendix A. 
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Hence we observe that, the MHV polygons may have a special reduction formula where all the 
coefficients of lower-order scalar loop integrals can be generated solely by the one that corresponds 
to its unique highest-order scalar loop integral. For this particular 5-pt case, this structure comes 
right from the unique reduction formula of MHV pentagons. 

To further appreciate the power of this Q functional coefficient, we go on to 1-loop 5-pt gravity 
amplitude. Since the structure of the expansion onto loop integrals satisfies the full color-kinematics 
duality, by the double-copy construction originally proposed by Bern, Carrasco and Johansson [ :], 
we can immediately substitute the color factor by another copy of the kinematic factor in each term 
in the the expansion, and the resulted formula is expected to be the correct gravity amplitude, which 
has already been conffimed in [8]. Moreover, it is not hard to check that in 5-pt case, the resulted 
expression can even be directly re-summed to be a formula purely consisted of MHV pentagons 
again 



5 From Yang-Mills to Gravity: 6-pt 1-Loop Amplitudes 



Since in 5-pt case we have observed that the MHV polygons serve as a very nice basis for A/" = 8 
supergravity amplitude, we would like to see whether they continue to work for more particles. But 
we need to start by checking 6-pt case, since redundancies start to occur in the reduction formula 
of MHV hexagons. The strategy is still to expand the MHV polygons into loop integrals, and 
seek for color-kinematics duality, and once this duality is satisfied, we can directly check whether 
double-copy construction gives the correct gravity amplitude. However, one needs to be cautious, 
because due to the redundancies the expansion does not in general have the correct color structure 
(in the sense that it matches with the corresponding trivalent diagram, upon which color-kinematics 
duality is based [i, 2, 8]. This is always true at tree level, but one needs to take loop propagators 
into consideration at loop levels. For more detailed discussion, please refer to Appendix A). In 6-pt 
case it turns out that, once we tune the expansion to have the correct color structure under our 
construction, color-kinematics duality is just a subsequent outcome, although a new phenomena 
would arise at the same time, which we will discuss in later parts of the paper. 

5.1 Analysis of the Redundancies 

In Section 2, we have already obtained a reduction formula (2.9) for MHV hexagons, which is 
manifestly cyclic symmetric. However, (2.9) is still not nice enough to work with. Instead, we 



(4.8) 




(1) 
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choose a particular point for the a parameter, and add deviations upon it 



e(123456) 

' ■ ' „(234561) = -M»||1MM _ A(123456). 

We can see the original constraint (2.10) on the a parameter is automatically satisfied. The virtue 
of analyzing around this particular point will be clear in the next subsection. But still A(234561) 
has no definition, and since a(234561) and a(123456) are assumed to be related by permutations, 
we would expect a(234561) can be naturally associated with a parameter A(234561) in the same 
way. And so it is necessary to further impose the condition 

(5.2) A(234561) = -A(123456). 

Then A{abcdef) shares the same symmetry as e{abcdef) under cyclic permutations. 

With these adjustments, the reduction formula of the MHV hexagons still only have pentagon 
integrals and box integrals, but no hexagon integral The hexagon integral should be introduced 
via the identities between hexagon integral and pentagon integrals, which can also be worked out 
systematically [2(i]. Since in the amplitude, each MHV hexagon is dressed by the color factor 
corresponding to the maximal color loop with the same sequence of labels, in a particular MHV 
hexagon we would add identity that involves only one hexagon integral whose labels fall into the 
correct sequence, otherwise in the resulted expansion formula the structure of the kinematic factor 
and the color factor would be drastically different, which is what we want to get rid of. Then take 
P( 123456) as an example, the identity that meets these requirements is unique 



(5.3) /6(1|2|3|4|5|6) - ^5(12|3|4|5|6) = 



We will dress it with a coefficient 7(123456). And by the consideration of the symmetry of the 
identity (5.3), this functional coefficient should satisfy 

(5.4) 7(123456) = 7(234561), 

(5.5) 7(123456) = 7(654321). 



5.2 Condition for Color-Kinematics Duality 

In order not to deviate into too much technical details, we choose to summarize here only the main 
results in the analysis of color-kinematics duality, and put the remaining details in Appendix B, C 
and D. 

^We don't consider here adding any tensor structures to the loop integrals, because by applying Carrasco and 
Johansson's ansatz [ ] with the most general tensor structures to 6-pt case, one can check that once the full color- 
kinematics duality is assumed to hold, any tensors of the type just vanish, and any tensors of the type P are 
purely gauge redundancies, which can be set to zero. We will not discuss this any further since it is not relevant for 
our current purpose. 
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As has been stated, in 6-pt case it is no longer true that the resulted formula automatically has 
the correct color structure. Since color-kinematics duality is based on the correct color structure, 
we should first make sure this condition is satisfied. We start by temporarily setting the parameter 
A = (so we only consider 7), and find that in the scalar loop integral expansion, box integrals of all 
types {l4{a\bc\d\ef), l4{ab\cd\e\f) and l4{abc\d\e\f)) already have the correct color structure. Since 
every MHV hexagon gives rise to its corresponding unique hexagon integral, the color structure of 
each hexagon integrals is also already correct. Then by looking at the coefficient in front of each 
pentagon integral (e.g. J5(12|3|4|5|6)), we obtain a constraint relating 7(123456) and 7(213456) 

7(123456) 7(213456) [12] [34] [56]si2S34g45g56(3|4 + 5|6)^ ^ 

^ ■ ^ (123456)e(123456) (213456)e(213456) (12)(34)(56)e(123456)e(213456)e(3456) 

As has been mentioned in the introduction, by permuting the labels, we can get constraints on 
other pairs of 7's related by transposition of two adjacent labels. So we can also regard (5.6) as a 
single functional constraint equation that depends on the choice of label ordering. 

Now we turn on the redundancies parameterized by A, and upon the previous analysis we only 
need to look at the effects of additional terms. Then the box integrals immediately constrain the 
A parameter to have the following form 

(5.7) Mahcdef) = — -^-j-niabcdef), 

e{abcdej) 

where n{abcdef) is completely symmetric under permutations of the labels, and so later on we 
will abbreviate it as k. Under this condition, the hexagon integrals are still untouched, while the 
constraint (5.6) from pentagon integrals now becomes 

7(123456) ^ 7(213456) 



(123456)e(123456) (213456)e(213456) 
, [12] [34] [56] (3|4 + 5|6)2 (3|4 + 5|6)2 

(12)(34)(56)e(123456)e(213456)e(3456) '^e(123456)e(213456)e(3456) ~ 

Since under (5.8) the correct color structure is guaranteed, we can go on to check the condition 
for color-kinematics duality. Interestingly, it turns out that the kinematic dual identities between the 
kinematic numerators of the hexagon integrals and pentagon integrals are exactly equivalent to (5.6) 
or (5.8) (depending on whether we set /t = 0), and the dual identities between kinematic numerators 
of the pentagon integrals and box integrals also automatically hold as a result of (repeatedly) 
applying (5.6) or (5.8). And the occurrence of k doesn't break the duality at all. In other words, 
starting from the Yang-Mills formula (3.5) in 6-pt case and introduce redundant parameters 7 
and K in the way as described above, just by requiring that the formula should have the correct 
color structure in its expansion onto scalar loop integrals, we know that these parameters only 
need to satisfy a single functional constraint (5.8) in addition to their own symmetries, and then 
color-kinematics duality just comes for free, almost the same as what happens at 5 points! 
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5.3 A Global Inconsistency 

Before going on, let us perform a consistency check on the constraint (5.6) or (5.8). We can see 
that, for every constraint with a specific configuration of the particle labels, it relates exactly two 
7's which differ by exchanging only one pair of neighboring labels. This can also be interpreted as, 
whenever we exchange two neighboring labels in 7, it would acquire some additional contribution as 
given by the inhomogeneous term in this constraint. Since this operation of exchanging labels allow 
us to exhaust the entire label configuration space starting from any specific point in it, we may 
expect that if we start from a certain point and travel step by step, and if in the end we travel back 
to the same starting point (label configuration), all the additional contributions acquired during 
the middle should add up to zero. This is an important consistency check in guaranteeing that the 
entire set of 7 parameters (with all in-equivalent label configurations) have a solution in terms of 
ordinary kinematics data. 

However, the constraint (5.6) or (5.8) that we have obtained in the previous subsection seems to 
be self-contradictory under this check. In more detail, the procedure as described above should be 
divided into two different categories. We may imagine the label arguments of 7 to reside on a circle, 
since by definition 7 should be invariant under cyclic permutation of the labels, and so there is no 
particular origin with respect to this operation. Then within the first category, no matter what 
kind of operation we do during the middle, after we have returned to the starting point, if the net 
effect is equivalent to that no labels have been brought to travel around the circle, we would find 
that all the additional contributions ultimately cancel out and the 7 just returns to itself, which is 
consistent. An example is illustrated as follows 

3 2 3 1 1 3 1 2 2 1 

4(^1 — > 4Q2-^ 4(^2 — > 4Q3^ 4Q3 

Q\ 56 56 56 56 56 

^ ' ^ 2 3 3 2 

- 4O1- 4O1. 

5 6 5 6 

However, in the other category, as long as the net effect is equivalent to that at least one label has 
been brought around the circle (i.e. let it to acquire a non-zero winding number), we would find a 
net discontinuity between the original 7 and the final 7 

(5.10) 1 {dbcdef ){ina_i = 7 (a&cde/) original + (additional term). 

The simplest example of this kind is 

32 31 13 43 43 43 

(5.11) 4Ql^ 4Q2— > 4(^2^ lQ2^ 5Q2^ 5Q2. 

56 56 56 56 16 61 

Attention should be paid that, no matter whether the additional term in (5.10) arise or not, its value 
depends only on the net effect (the winding numbers the labels acquired during the transformations), 
but not on the specific procedure of intermediate steps, and whenever a non-zero discontinuity is 
there, there is no means to tune k, in order to eliminate it. 
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This indicates that, if we regard the 7 and k as being constructed with ordinary variables (such 
as with the usual spin-helicity formalism) and attempt to solve the entire set of constraints, we 
would end up in having no solution for them. However, this inconsistency in 7 is not so harmful 
as it appears to be. The reason is that, back in the beginning when 7 is first introduced, it comes 
with a purely mathematical identity which is completely transparent to physical analysis, and 
characterizes an auxiliary parameter space; especially, this parameter contributes nothing to any 
leading singularities. If we do not insist that this parameter be constructed by ordinary variables, 
the inconsistency is not really there. Although not globally solvable, these "virtual" parameters 
do encode physical data in an astonishingly neat and compact way via the functional constraint 
on them, and allow us to travel between scattering amplitudes in different theories. Or put it in 
another way, the auxiliary space here only serve as a context where the constraint (the only object 
which carries physical information) can be properly described. These would be clarified in detail in 
the upcoming analysis. 



6 New Formula for 6-pt 1-Loop MHV Amplitude in A/^ = 8 
Supergravity 



From the previous section we know that although the existence of 7 and k in terms of ordinary 
kinematics data is at risk, the color-kinematics duality is very well satisfied, which still strongly 
tempts to suggest a "gravity" counterpart. Here in this section we are going to show that indeed 
we can still follow the original double-copy construction to obtain the correct formula for the 
corresponding 6-pt 1-loop gravity amplitude. This again confirms that these auxiliary parameters 
are not "locally" inconsistent. Instead, we should take the view that the constraints on these 
parameters work locally instead of globally, because all the physical information to be extracted 
only depends on certain subset of the constraints (or the same functional constraint with different 
arguments) which are mutually consistent with each other. 

The same as in 5-pt case, in 6-pt case the formula obtained by double-copy construction can still 
be reformulated purely in terms of MHV hexagons, which is equivalent to the corresponding Yang- 
Mills amplitudes with the color factor in front of each MHV hexagon substituted by the kinematic 
coefficient of the hexagon scalar loop integral with the same ordering of labels in the loop integral 
expansion 



In order to check the consistency of (6.1) with the correct quadruple cuts in gravity, we only need 
to check two types (since we are dealing with MHV polygons, leading singularities in quadruple cuts 
of the type {a\b\cd\ef) are trivially zero). Again we first assume k = 0. The first type is {a\hc\d\ef). 
Without loss of generality, we can assign particular configuration of the labels. So in this case we 



(6.1) 
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study the cut (1|23|4|56), which involves the following contributions 

7(123456) 7(123465) 7(132456) 7(132465) 
^'^ (123456)2 ^ (123465)2 ^ (132456)^ ^ (132465)^ ' 

The amazing fact is, we are still able to use the constraints (5.4), (5.5) and (5.6) repeatedly to fully 
eliminate the appearance of 7, and the final result is 



(6.3) 



[23]2[45][61]g23g45g56g6l(4|5 + 6|l)2 

(23)2(45)2(61)2(46)(51)e(234561)e(324561) 

[56][12][34]S56S12S23S34(1|2 + 3|4)2 [56][13][24]S56S13S32S24(1|3 + 2|4)2 



(123465) (56) (12) (34)e(561234)e(1234) (132465) (56) (13) (24)e(561324)e(1324) ' 



And this can actually be verified to be identical to the correct leading singularity of the gravity 
amplitude. 

For the other type we choose to study (1|2|3|456), which receives contributions from 6 terms 

7(123645) 7(123654) 7(123465) 7(123564) 7(123456) 7(123546) 
^ ■ '* (123645)2 ^ (123654)2 ^ (123465)2 ^ (123564)2 ^ (123456)2 ^ (123546)2 ' 

For this expression, we can use the same constraints in different ways to obtain different but equiv- 
alent final expressions, and still completely independent of 7, which can be shown to be equivalent. 
One example is 

/e(123456) _ [12]e(3456) \ [45][61][23]s45g6igi2g23(6|l + 2|3)2 
V (123456) (12) (34) (65) (36) (45) J (45)(61)(23)e(123456)e(123546)e(6123) 
[12]e(3456) [56] [12] [34]s56Si2S23^534(l|2 + 3|4)2 



(6.5) 



(12) (34) (65) (36) (45) (56) (12) (34)e(561234)e(651234)e(1234) 

[12]6(3546) [46][12][35]346gi2g23g35(l|2 + 3|5)2 

(12) (35) (64) (36) (54) (46) (12) (35)e(461235)e(641235)e(1235) 

[46][5l][23]g46g5lgl2g23(5|l + 2|3)2 [56][41][23]S56S41S12S23(4|1 + 2|3)2 



(123645) (46) (51) (23)e(123465)e(5123) (123654) (56) (41) (23)e(123564)e(4123) ' 

which again has been checked numerically to be equivalent to the correct leading singularity. 

Now we would also like to know how the k parameter would affect the leading singularities in 
the potential gravity amplitudes. The answer is that it has exactly zero effect. More explicitly, the 
quadruple cut (1|23|4|56)) receives an additional contribution 



(6.6) 



[ g56gl2g23g34(l|2 + 3|4)2 356^13^32^24(113 + 214)2 

{ (123465)e(561234)e(1234) (132465)e(561324)e(1324) 
[23]e(4561) 323345^56361 (4|5 + 6|1)2 I 

(23) (45) (61) (46) (51) e(234561)e(324561)e(4561) J 
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which can be checked to be exactly zero. And the cut (1|2|3|456) receives an additional contribution 

/ e(123456) _ [12]e(3456) \ ^45^61^12^236(6128) 
V (123456) (12) (34) (65) (36) (45) J e(123456)e(123546)e(6123) 

[12]e(3456) S56Si2S23S34(l|2 + 3|4)2 



(6.7) 



(12) (34) (65) (36) (45) e(561234)e(651234)e(1234) 

[12]6(3546) ^46^12^23^35(112 + 315)^ 

(12) (35) (64) (36) (54) e(461235)e(641235)e(1235) 

S46S5lSl2S23(5|l + 2|3)2 S56S4lSl2S23(4| 1 + 2 13)2 



(123645) (123465)e(5123) (123654)e(123564)e(4123) 

which also vanishes. As a result, in both types of quadruple cuts, the terms that contain k just 
cancel out completely. Hence we do not only get a formula (6.1) with the correct leading singularities 
of gravity amplitudes, but actually a one parameter family of them. The price for this is that we 
have to build the color-kinematics duality on a virtual auxiliary space. 



7 Discussions 

As a brief summary, we used MHV polygons to build 1-loop amplitudes for any number of particles 
in J\f = 4 super Yang-Mills theory, and up to 6 points in A/" = 8 supergravity, all of which are 
expressed as a simple sum over terms enjoying manifest cyclic and reflection invariance and are 
related purely by permutations. In other words, with this basis, the amplitudes can be generated 
with a single functional object (either a functional coefficient or a functional constraint). 

We have observed that non-trivial structures start to appear at 6 points. The coefficients 
multiplying the MHV hexagons cannot be expressed in terms of ordinary kinematics data. Instead, 
the amplitude here is formulated by imposing upon these coefficients a set of linear constraints, 
all of which have the same form (or, there is just a single functional constraint), which totally 
determine all the physical information needed from the amplitude. The striking fact is that, since 
these redundant parameters are not globally solvable, we haven't really obtained the real color- 
kinematics duality on the Yang-Mills side in our specific construction. However, the result shows 
that we can still circumvent it to get the correct gravity amplitude, just by imposing the algebraic 
relations but not solving them. By origin, 7 and k are only parameters that describe an auxiliary 
space. And so they can be regarded as virtual, and their only purpose is to form a context to 
formulate the constraint. It is really the constraint among these auxiliary variables that generates 
all needed physical data. 

As a natural generalization to the work presented in this paper, we can directly boost the form 
in (6.1) to a conjecture for any number of particles 

(7.1) M'i>= J: 

— (0"lO"2 • • • 
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with the 7(o"icr2 • • • <^n) still fully cyclic and reflection invariant but in general not solvable. Instead, 
we need to find out the constraints forced onto these coefficients, which are mutually consistent 
with respect to each quadruple cut but not necessarily consistent as a whole. By analogy, the 
coefficients here are still purely redundant parameters coming from adding kinematic identities 
into the expansion of MHV polygons in analyzing the color-kinematics duality (here by "purely 
redundant" we mean the parameters are dressed with pure zeros). We have seen that given the color- 
kinematics duality is satisfied, the gravity amplitude can be obtained just simply by substituting 
the color factor in front of each MHV polygons with the coefficient of the highest order scalar loop 
integrals in the expansion of the polygons. But in order to achieve that, at this point we still have to 
go back to the complicated expansions in terms of the scalar loop integrals and carefully fine-tune 
the redundancies. Especially, we have to seek for a manifestly symmetric decomposition of the 
MHV polygons into scalar loop integrals by introducing pure redundancies in a proper way, which 
in general seems hard. But most of the information in this analysis is actually not needed in the 
final answer. What we expect instead is a method that can work out the single coefficient needed 
to put in front of the MHV polygons or the constraints they shoTild satisfy. withoTit really working 
with the full set of loop integrals. Moreover, although the entire set of constraints are inconsistent, 
they have some good behaviors under the operations as described before, which seems to suggest 
there is some non-trivial mathematical structure to be better understood for this auxiliary space. 
And we hope that once this auxiliary space is described in a better way, it may help reach a general 
formula for 1-loop MHV amphtudes in J\f — S supergravity. 
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Appendices 

A A Brief Review of Color-Kinematics Duality 

In brief, color-kinematics duality refers to the situation that in a formula where the Yang-Mills 
amplitude is expressed totally in terms of trivalent diagrams (whose structures have one-to-one 
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correspondence to the color structures in terms of fabc), for every cofor Jacobi identity 
(A.l) Ci±Cj±Ck = 0, 

there exists a dual kinematic identity 

(A. 2) rii ±nj ±nk = 

sharing the same algebraic structure, but with each color factor substituted by its corresponding 
kinematic numerator [ ] (Notice that this is not necessarily a rational function. The remaining 
kinematic factor is a product of propagators indicated by the trivalent diagram). 

At tree level, when this duality is satisfied, the amplitude can in general be written in the form 



(A.3) = y: 



CiTti 



(lis) 



where the summation is over all different trivalent diagrams. As the simplest example, the 4-pt 
tree-level Yang-Mills amplitude can be written as (without the helicity factor) 

^-j ^(0) _ /l2a/a34 [12] ^ fnafa2i [13] 



si2 (23) (34) (41) (32) (24) (41)' 

which already satisfies the color-kinematics duality, with the dual kinematic identity associated to 
its unique color Jacobi identity 

fA5) [12] [13] _ 

^ ' (23) (34) (41) (32) (24) (41) 

4-pt is special in that, the duality holds for whatever form the amplitude is written in, which is due 
to the possibility of adding a redundancy parameterized by a 

(A.6) ^(°) = /i2a/a34 ( ^ ^2) (23)^) (41) ^ ") ^ ( (13) (32)^(24) (41) ~ ~ ^'^'^^''"^ 

while preserving the validity of the dual kinematic identity. For higher number of particles the 
formula for any amplitude does not automatically satisfy this duality, but the formula that meets 
this requirement is in general not unique. 

When a formula for tree-level Yang-Mills amplitude satisfies color-kinematics duality, Bern, 
Carrasco and Johansson conjectured that the correct formula for gravity with the same number of 
particles can be directly obtained by substituting in (A.3) each color factor Cj by another copy of 
its corresponding kinematic numerator rij [1] 



(A.7) = y: 



njUi 

WW- 



This is called double-copy construction. Especially, the two copies of kinematic numerator can come 
from two different formulas that satisfy color-kinematics duality respectively. This conjecture has 
been proved at tree level [■'>]. 
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The story is largely the same at loop level, except that here the duality is claimed to hold down 
to the integrand [_]. In other words, the amplitude can be expressed as 



(A.8) 



Cini{{p}, {k}) 



where the summation is still over all different trivalent diagrams, {k} is the set of loop momenta to 
be integrated over, Y[ is the product of loop propagators, and H is the product of propagators 
corresponding to additional tree structures of the trivalent diagram (if there is any). Note that, at 
loop level the kinematic numerator ni{{p}, {k}) is in general a function of both the external on-shell 
momenta and the loop momenta, so it may acquire tensor structures (in terms of loop momenta). 
The way to go to gravity is still to substitute each color factor Cj by another copy of the kinematic 
numerator nj({p}, {k}) 



(A.9) Ml'^ = J2l {Ud'k 



However, if it is known that any tensor structures appearing in rij can be gauged away (such as 
in 5-pt case as discussed by Carrasco and Johansson), then we can bring the color factor together 
with its corresponding kinematic numerator out of the integration in (A.8) 

where / denotes the ordinary scalar loop integral corresponding to the trivalent diagrams (without 
the propagators from the tree structures). And the conjectured gravity amplitude is 

Since in the cases that we study here the tensor structures are still redundant, we always use (A. 10) 
and (A. 11) when analyzing color-kinematics duality. 



B Analysis of the Color Structure 

This appendix summarizes the detailed analysis on the condition for 6-pt Yang-Mills amplitude to 
have the correct color structure when MHV hexagons are expanded into scalar loop integrals. To 
simplify analysis, we temporarily set the parameter A = 0. Then in the hexagon reduction formula, 
box integrals of the type l4{a\bc\d\ef) are 

^ ' ^ e(1234)e(4561) ^ ' ' ' ' 



17 



In order to obtain the correct color factor for this type of box integral, since this color factor is 
a linear combination of four different color factors of the highest level, we should expect that in 
the expansion of the full amplitude (multiply the above by the corresponding Parke- Taylor form 
and then sum over MHV hexagons), the kinematic factors in front of l4^{a\bc\d\ef), l4{a\bc\d\fe), 
l4^{a\cb\d\ef) and l4^{a\cd\d\fe) should be equivalent up to a possible minus sign, and this minus 
sign occurs when the two integrals differ by exchanging labels only within either (be) or (e/). It 
is straightforward to check that (B.l) meets this requirement, so this type of box integrals already 
have the correct color structure at the particular point as chosen in (5.1) when A = 0. And then it 
is also easy to see the same goes with the box integrals of the type l4^{ab\cd\e\f), which are 

(B.2) _^[12](23)|34](46)|56K61)...(5|3H-4|6)^ 

^ ' ^ e(3456)e(5612) ' ' ' ^ 

Then for the third type l^{abc\d\e\j\ the combined formula is 

^45^56(34) [46] (61) / [12](23)[34](45)[56][3|4 + 5|6) 
^ e(123456) V e(3456) 

^(12)123][45](56)[61](4|5 + 6|l|X 
e(4561) / 

The analysis is a little bit different for this case because of the color structure. For example, this 
particular box integral /4(123|4|5|6) would give contributions to two color structures faihfm-i and 
fa3bfbi2- So for each coefficient of this type of box integrals, we need to separate it into two parts 
and compare each part with the corresponding one from other box integrals that would contribute 
to the same box color structure, which is hard to be carried out. However, a careful look at the 
decomposition of color factors shows that another box integral 14(321 14|5|6) exactly and only gives 
contributions to the same two color structures, and with the same sign. By the fact that for each 
color structure, the decomposed coefficient should match between the two box integrals, the two 
entire coefficients should also be equivalent, and this has been verified by numerical check. Although 
this hasn't fully shown that the correct color structure can be obtained for this type of box integral, 
it is already a strong consistency check. And later in the analysis of color-kinematics duality, we'll 
further strengthen this argument. 

Then what is remaining are the pentagon integrals, which will give non-trivial constraints. The 
hexagon and pentagon integrals in the reduction formula are 

7(123456)/6(1|2|3|4|5|6) 



In order that the pentagon integrals have the correct color structure, we should have 
(B.5) 



r [12](23)[34](45)[56](61).34^45^56(3|4 + 5|6)^ -,Y123156) '^^^^^^ 1 ^ 
1 e(123456)e(3456) e(123456) J (123456) 



r [2 1] ( 1 3) [34] (45) [56] (62) .34^45 ^56 (3 1 4 + 5 1 6) ^ _ -, (213156) '^^^^^^ 1 ^ 
\ e(213456)e(3456) ^ e(213456) J (213456) ' 
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which gives 



,g . [12][34][56]gi2534g45g56(3|4 + 5|6)^ 7(123456) 7(213456) ^ 

^ ■ ' (12)(34)(56)e(123456)e(213456)e(3456) (123456)e(123456) (213456)e(213456) 

The constraints on 7 with the other sequences of labels are related just by permutations. So we 
can regard (B.6) either as a single functional constraint equation, or together with its permutations 
forming a set of constraints on the set of 7's with all possible orderings of particle labels. 



C Effects of the Parameter A 

This appendix complements Appendix B in studying the effect of A on the color structures. When 
A is non-zero, the terms corresponding to the type l4{a\bc\d\ef) receive additional contributions as 

Since the original part itself already gives the desired color factor, if this new expression should 
preserve this property, the additional term containing A should satisfy 

si4e(123456)(l|2 + 3|4)2 A(123456) _ Si4e(132456)(l|3 + 2|4)2 A(132456) _ ^ 
^ ' ^ e(1234)e(4561) (123456) e(1324)e(4561) (132456) ~ ' 

which is simplified to 

e(123456)A(123456) e(132456)A(132456) 



(C.3) 



(123456) (132456) 



This already gives a very strong constraint on A. But to further confirm it, take another type of 
box integral l4^{ab\cd\e\f) 

Again by requiring the desired color factor for this box integral, we would get 

e(123456)A(123456) e(213456)A(213456) 



(C.5) 



(123456) (213456) 



We may observe that given the symmetry of A, (C.5) is compatible with (C.3), and they together 
indicate that the particular form 

e{abcdef)A{abcdef) 
(abcdef) 
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should be invariant under all permutations of the labels. In terms of the new parameter k, terms 
involving the type Ii{abc\d\e\f) acquires the additional contribution 

(r7, ^ 545^56 [46] ( 123456) / (34) [3 14 + 5 1 6) (61)(4|5+^\ 

^^■'^ e(123456) [ e(3456) e(4561) ) ^^(123|4|5|6). 

Numerically, take J4(123|4|5|6) for example, its additional kinematic factor is invariant under the 
exchange of label 1 and 3, indicating that k doesn't break the correct color structure for this type 
of box integral as well. 

The pentagon integrals in the reduction formula now takes the form 



(C.8) ^« 



r[12](23)[34](45)[56](61)s34S45S56(3|4 + 5|6)= 



e(123456)e(3456) 

(123456).».«..a(3|4 + 5|6)'^ ^ .(3456) U^(,,|3|,|,|e) 

e(123456)e(3456) ^e(123456) J I I I I ^ 



Due to this additional term involving k, if we still want to preserve the correct color structure for 
pentagon integrals, the original constraint (B.6) on 7 should be modified to 

[iZ\ [64:\ [OOJS12S34S45S56 (3|4 + 5|6)2 S12S34S45S56 (3|4 + 5|6)2 

K- 



(12)(34)(56)e(123456)e(213456)e(3456) e(123456)e(213456)e(3456) 
^ ■ ^ 7(123456) 7(213456) 

^(123456)e(123456) ^ (213456)e(213456) ~ 



D Analysis of Color-Kinematics Duality 

This appendix studies the condition for color-kinematics duality when the scalar loop integral 
expansion of 6-pt Yang-Mills amplitude we have obtained has the correct color structure. 

Again let us set k = 0, and we first look into the relations between hexagon integrals and pen- 
tagon integrals. In order that color-kinematics duality holds at this level, for example, from the pen- 
tagon integral /5(12|3|4|5|6), together with two hexagon integrals /6(1|2|3|4|5|6) and /6(2|1|3|4|5|6), 
we should have the following constraint 

r [12] (23) [34] (45) [56] (61).34^45^56(3|4 + 5|6)^ _ ,^23156) '^^^^^^ ] ^ 
\ e(123456)e(3456) e(123456) J (123456) 

^ ' ^ _ ( 7(123456) 7(213456) 1 1 

~ \ (123456) (213456) J s^^' 

Then by using the identity 

(D.2) e(123456) + e(213456) + Si2e(3456) = 0, 
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it is easy to show that (D.l) is actually equivalent to the constraint (B.6) on 7 from the analysis of 
color structure. 

For the box integrals of the type l4{a\bc\d\ef), we pick up J4(l|23|4|56) as an example, and check 
the relation between this and /5(1|2|3|4|56) and /5(1|3|2|4|56). Substraction of the coefficients of 
the two pentagon integrals gives 

(D.3) 



/ [56](61)[12](23)[34](45)gi2g23g34(l|2 + 3|4)^ e(1234)7(561234) | 1 

\ e(561234)e(1234) e(561234) J (561234) 

/ [56] (61) [13] (32) [24] (45)^13^23^24(112 + 3|4)^ 6(1324)7(561324) ) 1 

\ e(561324)e(1324) e(561324) J (561324)' 



which, after imposing the constraint (B.6) on 7 from the pentagon color factors, results in an 
expression that is entirely independent of 7 

, [56][12][34]si2S23S34(l|2 + 3|4)2 ^ [56] [13] [24]si3S23S24(l|2 + 3|4)2 



(56)(12)(34)e(561234)e(1234) (56)(13)(24)e(561324)e(1234) 
^ ■ ' [23][45][61]s23S45S56S6i(4|5 + 6|l)2e(1234) 



(23) (45) (61)e(561234)e(561324)e(4561) 

If color-kinematics duality is to hold between this type of box integrals and the pentagon integrals, 
we should expect the above expression to be equivalent to the coefficient of the box times S23 

(12)[23](34)[41](45)[56](61)[14](1|2 + 3|4)^ 
^ ■ ' e(1234)e(4561) (123456) 

This has been verified with numerical values. Given this is satisfied, the other dual kinematics 
relations in the same category are automatic just by permuting the labels. The interesting fact here 
is that we don't even need to use the explicit form of the coefficients of the hexagon integrals. 

Then go on to study the type l4{ab\cd\e\f). For this we pick up /4(12|34|5|6) together with 
/5(1|2|34|5|6) and /5(2|1|34|5|6). Substraction of the coefficients of the two pentagon integrals gives 



(D.6) 



/ [34](45)[56](61)[12](23)556g6igi2(5|6 + 1|2)^ 6(5612) .g^^g^,. 

1 e(345612)e(5612) e(345612)^^ V (345612) 

/[34](45)[56](62)[21](13)s56S62Si2(5|6 + 2|l)2 _^^^^^q2i) 



\ e(345621)e(5621) e(345621) (345621)' 

Similarly, after imposing the constraints (B.6), this reduces to another 7-independent expression 

[34][56][12]556g6igi2(5|6 + 1|2)^ [34][56][12]556g623i2(5|6 + 2|1)^ 
(34)(56)(12)e(345612)e(5612) (34)(56)(12)e(345621)e(5612) 
^ ■ ' [12] [34] [56]gi2g34g45g56(3|4 + 5|6)^6(5612) 

(12)(34)(56)e(345612)e(345621)e(3456) ' 

which has been numerically checked to be equivalent to the coefficient of /4(12|34|5|6) times S12 

[12](23)[34](45)[56](61).56(5|3 + 4|6)^ ^ 
^ ■ ^ e(3456)e(5612) (123456) 
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This confirms that with the condition of correct color structure, color-kinematics duality also holds 
between pentagon integrals and the box integral of this second type. 

For the remaining type l4^{abc\d\e\f), we choose to study /4(123|4|5|6), and we would like to 
check whether the kinematic factor of this box integral with color structure /123456 is consistent 
with color-kinematics duality in relation to the pentagon integrals. As has be mentioned before, the 
color structure /123456 contributes to two box color structures, so we should study the dual kinematics 
relations for both structures, and only the sum (or substraction depending on the definition of the 
color structure) of the two results will equal the kinematic factor of J4(123|4|5|6). 

For the first color structure faibfb23, its related pentagon integrals are /5(1|23|4|5|6) and J5(23|l|4|5|6). 
The substraction of the two is 



(D.9) 



n23](34)[45](56)[61](12).45^56^6i(4|5 + 6|l)^ 6(4561) . ' 

\ e(234561)e(4561) e(234561)^^ V (234561) 

n23](31)[14](45)[56](62)5i4g45^56(l|4 + 5|6)^ _ Al^^(23U5Q)\ ^ 
\ e(231456)e(1456) e(231456)^^ V (231456)' 



Here we are still able to fully ehminate the 7, but should use two constraints. In the end, we would 
obtain the following expression 

[23][45][61]s45556S6i(4|5 + 6|l)2 , [23][U][5Q]suS4,s,eil\^ + 5\Q)^ 



(23)(45)(61)e(234561)e(4561) (23)(14)(56)e(231456)e(4561) 
^ ■ ^ [12][34][56]gi2g34g45g56(3|4 + 5|6)^e(4561) [13][45][62]gi3545g56g62(4|5 + 6|2)^e(4561) 
(12) (34) (56)e(123456)e(213456)e(3456) (13) (45) (62)e(134562)e(314562)e(4562) 

For the other color structure fasbfbu, we can obtain the formula just by doing the permutation 
{1 — ?■ 3, 2 — 7- 1, 3 — ?■ 2} to the above expression, which is 

[12][45][63]S45556S63(4|5 + 6|3)2 , [12] [34] [56]s34545S56(3|4 + 5|6)2 



(12)(45)(63)e(124563)e(4563) (12)(34)(56)e(123456)e(4563) 
^ ■ ' [31][24][56]33ig24g45g56(2|4 + 5|6)^6(4563) [32][45][61]332g45^56^6i(4|5 + 6| 1)^6(4563) 
(31) (24) (56)e(312456)e(132456)e(2456) (32) (45) (61)e(324561)e(234561)e(4561) 

Then, since the color structures under study are decomposed as 

^j-j fabdfblcfc23 falbfb2cfc3d falbfb3cfc2d fa2bfb3cfcld ~l~ fa3bfb2cfcld} 

fabdfb3cfcl2 = fa3bfblcfc2d ~ fa3bfb2cfcld ~ falbfb2cfc3d + fa2bfblcfc3d, 

it is (D.IO) — (D.ll) that should be equivalent to the kinematic factor of the box integral /4(123|4|5|6) 
multiplied by S123 = (pi + P2 + P3Y 



^45^56 [46] / [12](23)[34](45)[56](61)(34)[3|4 + 5|6) 
e(123456)V e(3456) 
^ ■ ^ (12)[23](34)[45](56)[61](61)(4|5 + 6|1] ^ 
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This equivalence has also been confirmed numerically. In doing this, we've not only verified the 
color-kinematics duality for this case, but also obtained the desired decomposition of the kinematic 
factor in (B.3) into (D.IO) and (D.ll) for l4^{abc\d\e\f) in order to give rise to the desired color 
structure for this type of box integrals, thus filling in the gap left by previous analysis. 

So far we see that when n = 0, the color-kinematics duality is automatically satisfied. Now we 
turn on k, and check whether this parameter would break the duality. Here we should be more 
careful, because non-vanishing n would modify the constraint on 7, and so in checking the duality, 
wherever 7 appears, we should count the additional term in (C.9). 

We still check the examples as we choose before, and only focus on the additional terms that 
would appear. At the level between hexagons and pentagons, if we add the k term from the 
coefficient of the pentagon into the L.H.S. of the original relation (D.l), we can observe that the 
modified relation is again just the condition of correct color structure for pentagons with k turned 
on. So the duality at this level is un-touched. 

At the level between pentagon integrals and box integrals, we first go back to the type involving 
l4{a\bc\d\ef) (e.g. /4(1|23|4|56)). Now the coefficient of the pentagon integrals J5(l|2|3|4|56) and 
/5(1|3|2|4|56) together with the constraint on 7 would add in (D.4) the following additional terms 

gl2g23g34(l|2 + 3|4)^ ^^^ 513^32524(113 + 2|4)^ ^^^ ^23^45^56^61 (4|5 + 6| 1)^1234) 

^ ■ ^ e(123456)e(1234) e(132456)e(1324) e(234561)e(324561)e(4561) 

while on the box integral side, its coefficient (times S23) contains the additional term 

mi^^ ^i4(l|2 + 3|4)^ 

^^•''^ -e(1234M4561)"-'^^- 

We have verified numerically that (D.14) and (D.15) are indeed identical. 

For the type involving l4{ab\cd\e\f) (e.g. /4(12|34|5|6)), the pentagon integrals J5(l|2|34|5|6) and 
75(21 1|34|5|6) together with the constraint give additional terms in (D.7) 

. 55656l5l2(5|6+l|2)^ ^^^ g5656252l(5|6 + 2|l)^ ^^^ gl2534545556(3|4 + 5|6)^e(5612) 

^ ■ ^ e(123456)e(5612) ^ e(213456)e(5621) ^ e(123456)e(213456)e(3456) 

And numerically, we have also checked that this is equivalent to the additional term given by the 
corresponding box integral (times S12) 

(D17) 556(5|3 + 4|6)^ ^^^ ^ 

^ ■ ^ e(3456)e(5612) 

For the remaining type involving l4^{abc\d\e\f) (e.g. /4(123|4|5|6)), we first analyze the rela- 
tion from J5(l|23|4|5p and /5(23|1|4|5|6), which together with the constraint adds into (D.IO) the 
following terms 

_ 54555656l(4|5 + 6|l)^ ^^ 5l4545556(l K + 5|6)^ 

e(234561)e(4561) e(231456)e(1456) 

^ ■ ^ 5l2534545556(3|4 + 5|6)^e(4561) ^^^ gl3345556562(4|5 + 6|2)^e(4561) ^^^ 

e(123456)e(213456)e(3456) ^ e(134562)e(314562)e(4562) ^' 
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Then in another relation (D.ll), we would get the following terms 



_ g45g56g63(4|5 + 6|3)^ ^^^ 534g45g56(3|4 + 5|6)^ 

e(124563)e(4563) e(123456)e(3456) 

^ ■ ' 331^24^45^56(214 + 5|6)^e(4563) 332g4535636l(4|5 + 6|1)^6(4563) 

e(312456)e(132456)e(2456) e(324561)e(234561)e(4561) ^' 

Similarly, in order for this type of dual kinematics relations to hold, we should expect that (D.18) — (D.19) 
matches with the additional terms from the corresponding box integral (times S123) 

m 9m ^45^56 [46] / (34)[3|4 + 5|6) (61)(4|5 + 6|1] >^ 

e(123456) e(3456) e(4561) ) ' 

And we have also verified this numerically. 



E Transformation to Explicit Formula for 6-pt 1-Loop MHV 
Gravity Amplitude in Terms of MHV Hexagons 

In this appendix we show that starting from our formula (6.1) for 6-pt 1-loop MHV amplitude in 
A/" = 8 supergravity, which is constructed with MHV hexagons together with the virtual coefficients 
7, it is also possible to obtain a formula, which is expressed explicitly in terms of ordinary kinematics 
data and meanwhile enjoys correct residues on all contours. 

In order to achieve this, pay attention that in general there are many non-planar identities 
relating different MHV polygons. Specifically for MHV hexagons, a direct analysis on leading 
singularities shows that, among the 60 inequivalent MHV hexagons, the maximally independent 
set only contains 40, so we would expect a great many identities. The smallest of such identities 
contains 8 MHV hexagons, which are related by exchanging labels within each pair, if we regard all 
the 6 labels as being separated into 3 adjacent pairs. One example is 

^6(123456) - 7^6(123465) - ^6(124356) + ^6(124365) 
^ ■ '* -^6(126543) + 7^6(125643) + ^6(126534) - 7^6(125634) = 0. 

Since these are purely mathematical relations, we are allowed to freely add them with any coefficient 
into the amplitude and meanwhile preserve all necessary physical data. 

Now go back to our formula for 6-pt gravity amplitude. The trick is to modify the coefficient of 
each term appearing in (6.1), say for 7(123456) 

7(123456) _ [16] [23] [45] 7(612345) [12] [34] [56] 7(123456) 

^ ■ (123456)2 ~ "(16) (23) (45) (612345)e(612345) ~ (12) (34) (56) (123456)e(123456) ' 

and the rest can be obtained just by permutations. A straightforward combinatoric counting gives 
60 inequivalent MHV hexagons, so after the modification we have doubled them to 120. Then notice 
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that now in front of each term, we always have a pre-factor of the pattern 



[ab][cd\[ef] 
{ab){cd){ef)- 

And in the 6-pt case there are altogether 15 of them. Then it is easy for us to check that the 120 
terms in the modified formula exactly fall into 15 groups, each involving 8 terms, and with the same 
pre-factor. If we denote 

j(abcdef) = 'yi^'bcdef) 

(ahcdef)e{abcdef) ' 

then for example, one of the group is 



[7(123456)^123456 + 7(123465)^123465 + 7(124 3 56)^124356 + 7(124365)P: 



124365 



(E.3) 

[12] [34] [56] 
(12) (34) (56) 

+7(125634)Pi25634 + 7(125 6 43)^125643 + 7(126534)^126534 + 7(126 5 43)Pi26543] • 

But again, all 7 are just simply related by the constraint (5.6) (since here we are only interested in 
obtaining an explicit example, for simplicity we just set k = 0). For this group we can bring the 
other 7 to 7(123456) with some additional terms purely of ordinary spinors, and when we collect 
all of the coefficient of the remaining 7(123456), we get 

[12] [34] [56] :,..oo,.g.r^ ^ V +V 

~ ^12) (34) (56) ^^ '^ L/-^123456 - /-'123465 " /•^124356 + /•^124365 

—'^125634 + "^125643 + "^126534 " "^126543] • 

However, the linear combination of the eight MHV hexagons in the big square bracket is just one of 
the standard identities of these kind of mathematical objects as shown in (E.l) that should always 
vanish. Since the same situation can be shown to occur in all other groups, to this point we have 
shown the 7 factors are completely eliminated from the entire formula, and what is left over is just 
the 1-loop gravity amplitude in the ordinary spin-helicity formalism, despite the fact that now the 
expression would look quite complicated and asymmetric. 
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